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Aim To analyze the haplotype of the Ezhava population of 
Kerala, south India, using 8 short tandem repeat (STR) loci 
on the Y chromosome and trace the paternal genetic lin- 
eage of the population. 

Methods Whole blood samples (n = 104) were collected 
from unrelated healthy men of the Ezhava population over 
a period of one year from October 2009. Genomic DNA 
was extracted by salting out method. All samples were 
genotyped for the 8 Y-STR loci by the AmpFiSTR Y-filer PCR 
Amplification Kit. The haplotype and allele frequencies 
were determined by direct counting and analyzed using 
Arlequin 3.1 software, and molecular variance was calculat- 
ed with the Y-chromosome haplotype reference database 
online analysis tool, www.yhrd.org. 

Results Among the 104 examined haplotypes, we found 
98 unique ones. The average gene diversity was 0.669, 
with the highest diversity of 0.9462 observed for the bi- 
allelic Y-STR marker DYS 385. The allele frequency among 
DYS loci varied between 0.0096 and 0.75. Out of the 104 
haplotypes, 1 0 were identical to the Jat Sikh population of 
Punjab, which is the greatest number among the Indian 
populations, and 4 to the Turkish population, which is the 
greatest number among the European populations. Ac- 
cording to the allele frequency of Y-STR, the Ezhavas were 
genetically more similar to the Europeans (60%) than to 
the East Asians (40%). 

Conclusion The vast majority of haplotypes were ob- 
served only once, reflecting the enormous genetic hetero- 
geneity of the Ezhavas. Based on the genotype, the Ezha- 
vas showed more resemblance to Jat Sikh population of 
Punjab and theTurkish populations than to the East Asians, 
hence indicating a paternal lineage of European origin. 
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Due to the geographical position of the Indian Peninsula 
between Africa, the Pacific, and west and east Eurasia, dif- 
ferent populations have moved through its territory. This 
is why ethnic Indian population shows enormous cultur- 
al, linguistic, and genetic diversity. Indian tribal and caste 
populations derive heritage largely from the Pleistocene 
southern and western Asians, receiving limited gene flow 
from external regions since the Holocene (1). Also, Indian 
castes have been found to be more closely related to the 
Central Asians than to the Indian tribal groups (2). 

The long seacoast of Kerala on the southern-most part of 
India has provided a gateway to India for many Asian, Eu- 
ropean, and Srilankan missionaries and traders. Non-tribal 
communities of Kerala, as shown by a human leukocyte 
antigen (HLA) analysis, were influenced by Dravidian, Indo- 
European, and East Asian gene pools (3). The Ezhava pop- 
ulation of Kerala, according to the allele frequency distri- 
bution, had features of European, Central Asians, and East 
Asian gene pools. Mitochondrial DNA studies also validated 
the presence of two distinct, eastern and western Eurasian- 
specific lineage groups in India, suggesting that there were 
at least two separate migration events to India (4). 

Due to the unique biology of the Y-chromosome, its ge- 
netic markers have been used in many forensic and evolu- 
tionary studies to determine patrilineal relationships with- 
in and between population groups. It has been suggested 
that, due to different distribution of region-specific allele 
frequencies, Y-short tandem repeats (STR) can be used to 
compare closely related populations (5,6). Previous genetic 
studies on the Ezhavas of south India failed to achieve a 
consensus on their paternal origin. In view of these diverse 
opinions based on HLA polymorphism and mitochondrial 
DNA analysis, this study aims to collect conclusive genetic 
data for a better understanding of the paternal origin of 
the Ezhavas. We present the haplotype analysis of the 8 Y- 
STR loci included in the European minimal haplotype set in 
104 men from the Ezhava population to explore their ge- 
netic relationships with the European and East Asian popu- 
lations. This is the first report on the Y-STR profile in Kerala 
population. 

MATERIALS AND METHODS 
Study samples 

Whole blood samples (n = 104) were collected from unre- 
lated healthy male individuals of the Ezhava population of 
Kerala over a period of one year from October 2009. The 



familial histories of the participants were recorded to ex- 
clude related individuals before sample collection. Blood 
samples were collected using standard procedures in eth- 
ylene diamine tetra aceticacid (EDTA)-coated tubes. The in- 
dividuals gave their written informed consent and the ethi- 
cal approval was received from the Ethics Committee of 
the institution. 

DNA analysis 

DNA was extracted from EDTA blood samples with the salt- 
ing out method (7). All samples were genotyped for the 
8 Y-STR loci (DYS19, DYS385, DYS389I, DYS389II, DYS390, 
DYS391, DYS392, and DYS393) by the AmpFiSTR Y-filer 
PCR Amplification Kit (Applied Biosystems, Foster City, CA, 
USA). The Y-STR data of Turkish, German, and North Indian 
Jat Sikh and other Indian populations were obtained from 
previous studies (8-10) and the Y-STR haplotype reference 
database (www.yhrd.org). 

STATISTICAL ANALYSIS 

Allelic and haplotype frequencies were estimated by direct 
counting. Gene diversities were calculated using Arlequin 
3.1 software (11), according to the formula G.D = N (1-IPi 2 ) 
/ N-1 , where N is the population size and Pi is the allele fre- 
quency of the ;-th haplotype. Furthermore, we used YHRD 
(www.yhrd.org) in order to determine the similarity of the 
Y-STR markers with the other publicly available population 
data. The Ezhava population was compared with other In- 
dian populations and with selected world populations in 
order to investigate the pattern of paternal contributions. 
A total of 1890 samples were analyzed: German popula- 
tion sample with 685 haplotypes, Turkish population sam- 
ple with 160 haplotypes, Ezhava population (our sample) 
with 104 haplotypes, Haryana Jat population sample with 
91 haplotypes, Punjab Jat Sikh population with 108 hap- 
lotypes, Andhra Pradesh, India (Brahmin) with 109 haplo- 
types, Himachal Pradesh, India (Saraswat Brahmin) with 
61 haplotypes, Jammu, India (Saraswat Brahmin) with 61 
haplotypes, Jharkhand, India (Munda) with 68 haplotypes, 
Jharkhand, India (Sakaldwipi Brahmin) with 65 haplotypes, 
Madhya Pradesh, India (Kanyakubja Brahmin) with 78 hap- 
lotypes, Punjab, India (Balmiki) with 62 haplotypes, Tamil 
Nadu, India (Iyengar) with 67 haplotypes, and Tamil Nadu, 
India (Kuruman) with 67 haplotypes. Analysis of molecular 
variance (AMOVA) was used to establish the total variance 
among and within groups. AMOVA was done using fixa- 
tion index analysis and with 10000 permutations. From 
the AMOVA calculation, we obtained pairwise differ- 
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TABLE 1 . Y-chromosomal short tandem repeat frequencies in Ezhava population of Kerala* 



Y-STR marker Y-STR marker 





DYS 


DYS 


DYS 


DYS 


DYS 


DYS 


DYS 


DYS 




DYS 


DYS 


DYS 


DYS 


DYS 


DYS 


DYS 


DYS 


ID 


19 


385 


389I 


38911 


390 


391 


392 


393 


ID 


19 


385 


389I 


38911 


390 


391 


392 


393 


E01 


15 


13,16 


12 


27 


24 


10 


11 


12 


E73 


15 


15,171 


13 


29 


22 


10 


10 


12 


E02 


14 


14,18 


13 


30 


25 


10 


11 


14 


E74 


15 


11,14 


13 


29 


25 


10 


11 


13 


E04 


15 


16,17 


14 


31 


22 


10 


11 


12 


E75 


16 


11,15 


14 


31 


24 


10 


11 


13 


E05 


14 


13,16 


13 


29 


22 


10 


14 


11 


E77 


15 


15,171 


13 


29 


22 


10 


11 


12 


E06 


14 


14,19 


13 


28 


23 


10 


10 


14 


E79 


14 


10,18 


13 


29 


23 


10 


10 


14 


E09 


13 


13,17 


13 


32 


26 


11 


13 


13 


E80 


14 


15,19 


14 


31 


23 


9 


11 


12 


E 11 


16 


11,14 


13 


30 


25 


11 


11 


13 


E81 


15 


11,13 


14 


32 


25 


10 


11 


13 


E 12 


15 


15,17 


14 


30 


22 


10 


11 


12 


E82 


16 


13,17.1 


13 


29 


23 


11 


11 


13 


E 13 


15 


15,17 


14 


30 


22 


10 


11 


12 


E83 


14 


11,19 


13 


29 


23 


10 


10 


14 


E 15 


13 


12,17 


13 


30 


26 


10 


13 


13 


E84 


16 


11,14 


13 


30 


25 


11 


11 


13 


E 16 


15 


11,19 


13 


31 


22 


9 


10 


14 


E86 


16 


11,14 


13 


31 


24 


11 


11 


13 


E 17 


15 


15,16 


13 


29 


22 


10 


11 


12 


E87 


15 


7,18 


13 


29 


22 


10 


14 


12 


E 18 


14 


13,15 


12 


28 


22 


11 


14 


11 


E88 


16 


10,14 


13 


28 


25 


10 


11 


13 


E 19 


15 


15,17 


14 


30 


22 


10 


11 


12 


E89 


14 


13,19 


12 


28 


22 


10 


14 


11 


E20 


14 


14,17 


14 


31 


24 


10 


11 


14 


E91 


14 


15,19 


13 


29 


23 


10 


11 


12 


E21 


15 


15,17 


13 


29 


22 


10 


11 


12 


E92 


15 


16,17 


14 


30 


22 


10 


11 


12 


E23 


14 


15,17 


13 


31 


24 


10 


11 


14 


E94 


14 


14,18 


13 


30 


22 


10 


14 


11 


E24 


14 


15,19 


13 


29 


23 


10 


11 


12 


E96 


17 


11,14 


14 


32 


26 


10 


11 


13 


E25 


15 


12,14 


14 


32 


25 


10 


11 


13 


E99 


14 


13,16 


12 


29 


22 


10 


14 


11 


E26 


15 


12,14 


14 


32 


25 


10 


11 


13 


E100 


15 


15,19 


14 


30 


22 


10 


12 


13 


E27 


17 


13,19 


14 


31 


23 


10 


10 


14 


E101 


15 


11,15 


14 


32 


24 


10 


11 


13 


E28 


14 


11,14 


14 


31 


25 


10 


11 


13 


E102 


14 


14,18 


13 


29 


22 


10 


14 


11 


E29 


17 


14,15 


13 


29 


21 


10 


11 


13 


E103 


14 


14,18 


13 


31 


25 


10 


11 


14 


E32 


17 


11,14 


13 


30 


25 


11 


11 


13 


E104 


14 


14,18 


13 


30 


24 


10 


11 


13 


E33 


15 


15,16 


13 


29 


22 


10 


11 


12 


E105 


14 


13,17 


14 


31 


22 


10 


14 


12 


E36 


15 


11,13 


14 


32 


25 


10 


11 


13 


E106 


15 


13,17 


12 


28 


24 


11 


11 


12 


E37 


16 


11,12 


13 


30 


24 


11 


11 


13 


E107 


15 


11,14 


14 


32 


25 


10 


11 


13 


E39 


16 


11,14 


13 


31 


25 


10 


11 


13 


E108 


15 


15,17 


13 


29 


22 


10 


11 


12 


E41 


14 


13,17 


12 


28 


22 


11 


14 


11 


E110 


14 


15,17 


12 


28 


22 


10 


14 


11 


E42 


15 


11,14 


13 


30 


25 


10 


11 


13 


R1 


14 


15,17 


13 


32 


24 


14 


11 


15 


E43 


14 


11,19 


13 


29 


23 


10 


10 


14 


R3 


14 


15,17 


13 


32 


25 


10 


11 


15 


E44 


15 


12,14 


14 


32 


25 


10 


11 


13 


R4 


14 


15,17 


13 


0 


24 


10 


11 


15 


E46 


15 


15,17 


14 


30 


22 


10 


11 


12 


R5 


14 


15,17 


13 


32 


24 


10 


11 


15 


E47 


14 


14,18 


12 


28 


22 


10 


14 


11 


R7 


16 


11,14 


13 


30 


25 


10 


11 


13 


E49 


17 


11,14 


14 


30 


25 


10 


11 


13 


X 


15 


11,14 


13 


31 


25 


10 


11 


13 


E50 


14 


13,18 


12 


28 


22 


11 


14 


11 


E 14 


15 


11,14 


13 


29 


25 


11 


11 


12,13 


E51 


14 


13,17 


14 


31 


22 


10 


14 


12 


E38 


15 


11,13 


14 


31,32 


25 


10 


11 


13 


E52 


15 


11,14 


13 


30 


25 


10 


11 


13 


E53 


15 


11,14 


12,13 


29,30 


25 


10 


11 


13 


E54 


15 


13,16 


12 


29 


23 


10 


14 


14 


E60 


15,16 


11,14 


14 


31 


25 


10 


11 


13 


E55 


15 


12,14 


14 


32 


25 


10 


11 


11 


E22 


15 


16 


12 


27 


25 


11 


13 


13 


E56 


14 


11,14 


14 


31 


25 


10 


11 


13 


E31 


15 


14 


12 


29 


24 


11 


13 


12 


E57 


16 


11,12 


13 


30 


24 


11 


11 


13 


E35 


15 


15 


13 


28 


22 


11 


11 


12 


E59 


15 


13,16 


12 


27 


24 


10 


11 


12 


E40 


14 


12 


13 


29 


23 


10 


10 


14 


E61 


16 


11,12 


13 


30 


24 


11 


11 


13 


E48 


15 


14 


12 


29 


24 


11 


13 


12 


E62 


14 


14,18 


12 


28 


22 


10 


14 


11 


E63 


15 


14 


12 


29 


24 


11 


13 


12 


E66 


15 


11,14 


13 


29 


25 


10 


11 


13 


E64 


14 


12 


13 


29 


23 


10 


10 


14 


E67 


14 


14,17 


14 


32 


25 


10 


11 


14 


E65 


15 


15 


13 


29 


22 


10 




12 


E68 


15 


12,15 


12 


29 


24 


11 


11 


13 


E76 


16 


11 


13 


30 


26 


11 




13 


E69 


15 


14,19 


13 


30 


22 


10 


11 


12 


E78 


15 


16 


13 


28 


22 


10 




12 


E70 


14 


15,19 


13 


29 


23 


10 


11 


12 


E90 


15 


14 


14 


30 


23 


10 




12 


E71 


16 


11,15 


13 


31 


26 


10 


11 


13 


E95 


15 


18 


13 


30 


22 


10 




12 


E72 


15 


15,16 


13 


29 


22 


10 


11 


12 


E109 


15 


14 


14 


30 


22 


10 




12 



*E,R,X - Ezhava samples. 
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ence and genetic distance between populations (multidi- 
mensional scaling analysis). 

RESULTS 

The haplotypes of 104 unrelated men from the Ezhava 
population of Kerala were investigated for 8 markers in- 
cluded in the European minimal haplotype set (Table 1). 
The majority of haplotypes were unique (98/104). While a 
larger number of individuals showed monoallelic condi- 
tion at the DYS19, DYS389I, DYS389II, and DYS393 loci, a 
few samples (4/104) were unique, showing biallelic con- 
dition or duplication at the loci. A total of 12.5% (13/104) 



of the individuals showed monoallelic condition at the 
normally biallelic DYS385 locus. The variant allele DYS 385 
*1 7.1 was observed in 3 samples (Table 1). 

Allelefrequencies at DYS loci varied between 0.75 to 0.0096 
(Table 2). Gene diversity ranged from 0.3875 to 0.9462 (av- 
erage, 0.6671) and the most polymorphic single locus 
marker was DYS389II (Table 2). 

The haplotypes of the Ezhavas were compared with the 
haplotypes of other Indian populations (Table 3). The pair- 
wise difference analysis results were 0.0528 for Punjab Jat, 
0.0691 for the Turkish population, 0.1024 for the Harya- 



TABLE 2. Allele frequencies and gene diversity values at 8 DYS short tandem repeat loci of the Ezhava population of Kerala 



DYS Allele 3891 390 38911 


19 


393 


391 


392 


Genotype 


385* 


9 






0.0192 




7,18 


0.0096 


10 






0.75 


0.0865 


10,14 


0.0096 


11 




0.1154 


0.2308 


0.6827 


10,18 


0.0096 


12 0.1923 




0.3462 




0.0096 


11,12 


0.0481 


13 0.5288 


0.0192 


0.4038 




0.0673 


11,13 


0.0288 


14 0.2881 


0.3077 


0.1346 




0.1442 


11,14 


0.1731 


15 


0.4904 








11,15 


0.0288 


16 


0.1442 








11,19 


0.0288 


17 


0.0481 








12,14 


0.0385 


21 0.0096 










12,15 


0.0096 


22 0.3365 










12,17 


0.0096 


23 0.1346 










13,15 


0.0096 


24 0.2019 










13,16 


0.0577 


25 0.2692 










13,17 


0.0481 


26 0.0481 










13,17.1 


0.0096 


27 0.0385 










13,18 


0.0096 


28 0.1154 










13,19 


0.0192 


29 0.3077 










14,15 


0.0096 


30 0.2596 










14,17 


0.0192 


31 0.1538 










14,18 


0.0673 


32 0.125 










14,19 


0.0192 












15,16 


0.0288 












15,17 


0.0769 












15,17.1 


0.0192 












15,19 


0.0481 












16,17 


0.0192 


11 0.0096 












12 


0.0192 












14 


0.0577 












15 


0.0192 


16 0.0192 












18 


0.0096 


GDV f 0.6063 0.7602 0.7908 


0.6475 


0.6923 


0.3875 


0.5059 




0.9462 


The allele frequencies for DYS385 were calculated for the combination of two alleles. 



tGene diversity value. 
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na Jat, and 0.1785 for Germans. These values show that 
Ezhava population is distant from other Indian popula- 
tions and close to the Punjab Jat and Turkish population, 
as well as the Jammu population (Saraswat Brahmin), with 
a pairwise distance value of 0.0561, and the Jharkhand 
population (Sakaldwipi Brahmin), with a value of 0.0617 
(Figure 1). 

Similarities in the Y-STR data were found between the 
Ezhavas, north Indian, and Turkish population (Table 4). 
Out of the 104 haplotypes, 10 haplotypes were identical 
to Jat Sikhs, 7 to the German population, and 4 to the Turk- 
ish population. One particular haplotype was found to be 
common to all 4 populations, with a very high incidence 
(8/108) in the Jat Sikh population of north India. 

Allelic distribution of the Y-STR markers of the Ezhavas was 
compared with the European and East Asian populations 
(Figure 2). The histogram showed that 60% of the mark- 
ers had similar allele frequencies to the European popula- 

TABLE 3. Analysis of molecular variance pair-wise distances 
based on the fixation index (F ) values between the Kerala 
Ezhava population and other Indian populations 

Population samples* Kerala Ezhava 

Jammu (Saraswat Brahmin) 0.0573 

Jharkhand (Sakaldwipi Brahmin) 0.0653 

Tamil Nadu (Iyengar) 0.0671 

Himachal Pradesh (Saraswat Brahmin) 0.0693 

Punjab (Balmiki) 0.0713 

Tamil Nadu (Kuruman) 0.0734 

Andhra Pradesh (Brahmin) 0.0824 

Madhya Pradesh, (Kanyakubja Brahmin) 0.0840 

Jharkhand (Munda) 0.1321 
*Names in parenthesis refer to the YHRD designations (available from 
http://www.yhrd.org). 



tion, while only 40% showed similarity to the East Asian 
populations. The Y-STR alleles showed marked similarities 
to European populations when compared using European 
minimal haplotype. 

DISCUSSION 

This study found that the Ezhavas showed more genotypic 
resemblance to the Jat Sikh population of Punjab and the 
Turkish populations than to East Asians, hence indicating a 
paternal lineage of European origin. 

Figure 1. 



Haryana, India [Jat] 



Germany 



Punjab, India [Jat Sikh] 



Turkey 



0 Kerala, India 
[Ezhava] 



Dimension 1 
stress = 0.0273 



Multi Dimensional Scaling plot showing the relationship between Kerala 
Ezhava population with other Indian and European populations. 



TABLE 4. Comparison of haplotypes of Y chromosomal short tandem repeats in populations from Germany (G), north India (Nl), 
Turkey (T), and Ezhava population (KE) of South India 

Y-STR marker 



DYS 19 


DYS 385 


DYS 389 I 


DYS 389 II 


DYS 390 


DYS 391 


DYS 392 


DYS 393 


G 


Nl 


T 


KE 


15 


15,16 


13 


29 


22 


10 


11 


12 




1 




3 


14 


14,19 


13 


28 


23 


10 


10 


14 




1 




1 


16 


11,14 


13 


30 


25 


11 


11 


13 


1 


8 


2 


2 


15 


15,17 


14 


30 


22 


10 


11 


12 






1 


1 


14 


14,18 


12 


28 


22 


10 


14 


11 




1 




2 


15 


11,14 


13 


29 


25 


10 


11 


13 


1 






2 


15 


15,17 


14 


30 


22 


10 


11 


12 






1 


1 


15 


11,14 


13 


30 


25 


10 


11 


13 


2 






2 


16 


11,14 


13 


30 


25 


10 


11 


13 


2 


3 




1 


15 


11,14 


13 


31 


25 


10 


11 


13 




1 




1 
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Figure 2. 




19 389 I 389 II 390 391 392 393 



DYS markers 

Comparison of the DYS allele with the highest allele frequency between 
Kerala Ezhava - dark gray bars, European - closed bars, and East Asian 
population - light gray. 

To our knowledge, there are no published data about the 
genetic structure of Y-chromosome in Kerala population, 
although data on various other ethnic groups have been 
reported in the last decade (12-19). Y-STR haplotypes are 
useful for investigating and reconstructing the phylogeny 
of the more recently diverged paternal lineages (5,6), as 
well as for forensic/paternity testing (20-24). Based on HLA 
polymorphism analysis, Thomas et al (3) reported that the 
Ezhava population had features of European populations 
such as Belgian, German, and Scottish and that non-tribal 
communities of Kerala showed influences of Dravidian, In- 
do-European, and East Asian gene pools. Our study using 
Y-STR markers also confirmed their findings that the Ezha- 
vas derived their paternal lineage from European and East 
Asian populations. There are also similar reports that placed 
Indians closer to European populations than to either East 
Asians or Africans in the genetic distance trees (25). 

Duplications, though rare in the Ezhava population, have 
been observed in monoallelic markers like DYS1 9, DYS389I, 
DYS389II, and DYS393 loci in other populations (8,20,24). 
Duplications at DYS19 locus have been reported by previ- 
ous studies in the Turkish populations (8). Except one re- 
port (26) of the variant allele DYS385 1 7.1 in Croatian pop- 
ulation, the Ezhava population seem to be the only other 
population with variant DYS 385 1 7.1 allele. 

The present study established the European paternal lin- 
eage of the Ezhavas, as well as contributed to the creation 
of a database of 8 Y-STR markers in the population. A limi- 



tation of the study is the non-availability of comparative 
Y-STR-based data for other Kerala communities. Also, the 
study was limited to a particular geographic area with the 
predominant Ezhava community. The sample size was lim- 
ited to 104 participants as it was a pioneering study to es- 
tablish the sample size for a study with a larger sample 
size and wider geographic coverage. Future studies are 
planned to create a database for 17 Y-STR markers in all 
the communities in Kerala population and to enhance the 
sample number and geographic area of the study. 
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